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Abstract 

Context. Long soft gamma ray bursts (LGRBs) are usually associated with the death of the most massive stars. A large amount of 

core angular momentum in the phases preceding the explosion is required to form LGRBs. A very high initial rotational velocity can 

provide this angular momentum. Such a velocity strongly influences the way the star evolves: it is chemically homogeneously mixed 

and evolves directly towards the blue part of the HR diagram from the main sequence. 

Aims. We have shown that chemically homogeneous evolution takes place in the SMC, at low metallicity. We want to see if there is a 

metallicity threshold above which such an evolution does not exist. 

Methods. We perform a spectroscopic analysis of H-rich early-type WN stars in the LMC and the Galaxy. We use the code CMFGEN 

to determine the fundamental properties (T^s, L) and the surface composition of the target stars. We then place the stars in the HR 

diagram and determine their evolution. 

Results. We show that both the LMC and Galactic WNh stars we selected cannot be explained by standard stellar evolution. They are 

located on the left of the main sequence but show surface abundances typical of CN equilibrium. In addition, they still contain a large 

amount of hydrogen. They are thus core-H burning objects. Their properties are consistent with chemically homogeneous evolution. 

We determine the metallicity of the Galactic stars from their position and Galactic metallicity gradients, and conclude that they have 

0.6<Z<1.0. A moderate coupling between the core and the envelope is required to explain that stellar winds do not extract to much 

angular momentum to prevent a blueward evolution. 

Conclusions. We have shown that chemically homogeneous evolution takes place in environments with metallicity up to solar. In 

view of the findings that some long gamma ray bursts appear in (super-)solar environments, such an evolution may be a viable way to 

form them over a wide range of metallicities. 
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^D 1. Introduction they also lose a significant amount of angular momentum dur- 

Cn ing their life. But massive stars winds also depend on metallicity 

^H , Long soft Gamma-Ray Bursts (LGRB) are classically associ- (^^^ 2^).6-o.i iMokiemet all 120071 IVinket alJ 120011) : the lower 

^ , ated with the death of massive stars. The most convincing ev- j^e metal content, the smaller the radiative acceleration through 

idence for this is the clear association of some LGRBs with type Yxn&& and thus the lower the mass loss rate. The observation of 

^ . Ib/c supernovae (Galamaet al 1998; H joithet al. 2003 ). They i^^ger r otational velociti es in the LMC/SMC comp ared to the 

H ■ are also observed in star forming galaxies, often directly in sites G^i^xy dMartavan et al.1 120061: iMokiem et al.1 120061) is an indi- 

^: where stars are currently being formed (but see [Hammer et al. | ^ect confirmation of this effect (although the initial distribution 

|2006|). The conditions required to produce LGRBs impose that ^f rotational velocities might be different too). Massive stars in 

its progenitor must have an envelope free of hydrogen so the jet j^^.z environments are thus more likely to become LGRBs in 

can travel through it without being damped. The progenitor core t^e scenario wher e the angular momentu m evolution is governed 

must also retain a large angular momentum m order to produce ^y ^^^ss loss Isee lPetrovic et"d] (l2()05bl) for alternative possibil- 

a disk-jet structure around the black hole resulting from thecore jties).. But an unusually large initial rotational velocity is also 

collapse. The collapsar model proposed by W oosley et al. |(11993|) required since not all low-Z massive stars retain enough angular 

fulfills these conditions and is presently the main scenario in- momentum at the end of their Ufe to become a LGRB. 

voked to explain the existence of LGRBs. 

The conditions required by the collapsar scenario, especially Rapid rotation affects stellar evolution. iMaedeJ (1 19871) (see 

the large angular momentum of the core before explosion, have also iLanged Il992h showed that for very fast rotating massive 

implications for the properties of the progenitor star Since mas- stars, the mixing timescale becomes shorter than the nuclear 

sive stars lose mass at a l arge rate due to their strong radia - timescale. In that case, the material produced in the stellar core is 
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tively driven stellar winds (ICastor et al.lll975l iPuls et alll2008h . immediately mixed with the outer layers. The star is thus (quasi) 
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chemically homogeneous. The larger fraction of helium in the 
outer layers reduces the opacity, causing the star to become hot- 
ter. Consequently, chemically homogeneous stars evolve directly 
to the blue part of the HR diagram from the main sequence. 

We showed that quasi-chemically homogeneous evolution 
(CHE) was most hk ely followed by the star WR 1 in the SMC 
JMartinsetal.1 120091) . The analysis of its UV-optical-IR spec- 
trum indicated a hot temperature (65000 K), placing the star 
on the blue side of the main sequence in the HR diagram. At 
the same time, we determined that WR 1 still contained a large 
amount of hydrogen in its envelope. This was consistent with its 
spectral classification (WN5h) typical of WN stars showing the 
presence of hydrogen lines in their atmosphere. Classical evolu- 
tionary tracks are not able to reproduce both the position in the 
HR diagram and the large hydrogen content. When they reach 
the position of WR 1, they have long burnt all of their hydro- 
gen. Only chemically homogeneous evolution can explain the 
properties of WR 1 . This is strengthened by our determination 
of the CNO content, fully consistent with the values expected 
from CNO equilibrium. SMC WR 1 is thus very likely a massive 
star in the core-H burning phase and evolving homogeneously 
towards the blue part of the HR diagram. It validates the con- 
cept that some stars follow CHE even if this does not necessarily 
mean that they will end as a LGRB. 

This study was focused on stars in the SMC. One may won- 
der whether CHE al so happens at higher metallicity (see also 
iGrafener et al.l I2OIII) . This is important in the context of the 
collapsar model. The first determinations of the metallicity of 
LGRB hosts tended to show that low-Z environments were pre- 
ferred. For instance. iModiaz et al.l ( l2008b showed that type Ic su- 
pemovae associated with LGRBs were located in galaxies with 
metallicities systematically lower than type Ic SN without GRB. 
Various analysis of LGRB hosts con firmed that low m etallicity 
was preferred (Thone et al. 2008; Le vesque et al.ll201ll) . On the 
theoretical side. lYoon et al.l (120061) showed that stellar evolution 
at metallicities above ~0.6 ended with too low angular momen- 
tum in the core to produce LGRBs. A threshold in metallicity 
around 0.6 Z0 for the formation of LGRBs seemed to exist. 

Things have changed recently on the observational side with 
the discovery of at least two LGRBs at su per-solar metallicity 
dGraham et al.ll2009l: iLevesque etai]|20I0h . In addition, the de- 
termination of the mass-metalhcity relation for LGRB hosts in- 
dicates that solar and super-solar metallicities are not excluded, 
although the relation is shifted towards lower Z compared to the 
classical M-Z relation dMannucci et alj|201 ih . The tendency of 
LGRBs to be located in actively star-forming galaxies explains 
this offset: such galaxies are currently mainly low-mass galaxies. 

The question of the occurrence of CHE at high metallic- 
ity then arises. If massive stars never evolve homogeneously at 
high Z, this might favour alternative evol utionary paths to form 
LGRBs i n high metallicity env ironments. iPetrovic et al1(l2005bl) 
(see also IGrafener et all |20 1 2h showed that under the assump- 
tion of a weak coupling between the core and the envelope, the 
accreting star of a close binary, or a single star, might retain 
enough angular momentum after the red supergiant phase to pro- 
duce a LGRB. In the present paper, we report on the analysis of 
hydrogen-rich WN stars in the LMC and the Galaxy. In a man- 
ner similar to the WNh stars in the SMC, we determine their 
physical properties and conclude about their evolution. We show 
that CHE is likely to exist in the LMC and the Galaxy, at solar 
metallicity. 

The paper is organized as follows; in Sect. |2] we describe 
our sample and the observations; we then explain how we deter- 
mined the fundamental properties of the sample stars in Sect.|3j 



we then present our results in Sect. Hand discuss them in Sect. 
|5] The conclusions are gathered in Sect. |6] 

2. Sample stars and observations 

The optical data have been obtained from a variety of tele- 
scopes/instruments. Data f or the LMC stars hav e been collected 
at SAAO and are those of iFoellmi et al.l (l2003l). Det ails on the 
data reduction can be found in iFoelkni et al.l (l2003l) . We have 
used the average of all optical spectra in our analysis, to pro- 
vide average stellar parameters. The optical spectrum of WR 128 
was collected by one of us (LM) during run 089.D-0730(A)) on 
the 2.2m telescope equipped with FEROS spectrograph at the 
La Silla ESO observat ory. The reduction procedure is identi- 
cal to that described in ISana et all (l2006l) and was done under 
the MIDAS environment. Data for WR 7, WR 18 and WR 128 
were obtained with the spectrograph UVES on ESO/VLT (run 
080.D-0137(A)). UVES spectra were reduced using the UVES 
pipeline provided by ESO, which performs bias and inter-order 
background subtraction (object and flat-field), optimal extraction 
of the object (above sky, rejecting cosmic ray hits), division by 
a flat-field frame extracted with the same weighted profile as 
the object, wavelength calibration and rebinning to a constant 
wavelength and step, and merging of all overlapping orders. In 
addition to the optical spectra, we have collected UV data from 
the lUE archive. Information on the spectra are given in Table[T] 
When several lUE spectra are listed, we used the average spec- 
trum to perform our spectroscopic analysis. 

3. Modelling 

We have used the code CMFGEN jHillier & Milled [19981) to 
determine the stellar and wind parameters of the sample stars. 
CMFGEN solves the radiative transfer in the comoving frame 
in non-LTE conditions. A spherical geometry is adopted to cor- 
rectly take the wind extension into account. The rate equations 
and the radiative equation are solved simultaneously through an 
iterative scheme until convergence of the level populations has 
been achieved. The temperature structure is set by the radiative 
equilibrium constraint. The velocity law (and consequently the 
density law through the mass conservation equation) is the re- 
sult of the combination of an inner structure connected to a /? 
velocity law in the outer atmosphere. The inner structure is it- 
erated: once the level populations have converged, the total ra- 
diative acceleration is calculated and the momentum equation is 
solved to provide a new velocity law. This solution is connected 
to the y6 velocity law and the resulting velocity structure is used 
to compute a new atmosphere model with new level populations. 
A couple of these hydrodynamical simulations is performed. In 
the model computation, a constant microturbulent velocity of 50 
km s ' is adopted. The elements included in our models are H, 
He, C, N, O, Ne P, Si , S, Fe, Ni. We use the solar abundances of 
iGrevesse et al.l (|2007|) as reference. 

Once the atmosphere model is completed, a formal solution 
of the radiative transfer equation is performed to yield the emer- 
gent spectrum. A depth variable microturbulence velocity, from 
10 km s"' at the bottom of the wind to 10% of the terminal ve- 
locity in the outter layers is adopted. The emergent spectrum is 
subsequently compared to the observed data to determine the 
stellar and wind parameters. In practice, we proceed as follows: 

• Effective temperature: we use the ionization balance method 
to constrain T^s- The ratio of N iv to N v lines is a first indi- 
cator. N IV 4058, N IV 5205, N iv 5736, N iv 7 103, N v 45 19- 
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Table 1. Sample stars and observational information 



Star 


Optical spectrum 


Date of optical observation 


Resolution 


UV spectrum 


Date of UV observation 


Resolution 


Galaxy 


WR7 


VLT/UVES 


26-27 Jan 2008 


30000 


lUE sp03541 


6 Dec 1978 


250 


WRIO 


VLT/UVES 


26-27 Jan 2008 


30000 


lUE spl0748/spl3912 


4 Dec 1980/ 6 May 1981 


250 


WR18 


VLTAJVES 


26-27 Jan 2008 


30000 


lUE sp0933 


18 May 1980 


250 


WR128 


NTT/FEROS 


3 May 2012 


48000 


IUEspl5101 


26 Sep 1981 


18000 


LMC 


Bat 18 


SAAO/GS' 


28 Dec 1999 

16-20 Jan 2001 

1-7 Jan 2002 


900 


lUE sp38236/sp38247 


24 Fev 1990 


250 


Bat 63 


SAAO/GS' 


28 Dec 1999 

16-20 Jan 2001 

1-7 Jan 2002 


900 


lUE sp09157/sp09158 


31 May 1980 


250 




SAAO/GS stands for SAAO observatory / Grating spectrograph. Fron 








Notes. 1. 


ilFoellmiet all (120031) 





4523, Nv 4944 are the main indicators. We also use the 
strength of the He ii lines and, when observed, of He i 5876 
as a secondary temperature indicators. The Fev 1350-1400 
forest in the UV range provides an additional constraint. For 
Wolf-Rayet stars, a meaningful quantity is the temperature at 
an optical depth of about 20 (T*). It is appropriate for com- 
parison with temperatures provided by evolutionary models 
since it corresponds to the quasi-hydrostatic part of the at- 
mosphere. 

Luminosity: the luminosity is derived from the fit of the spec- 
tral energy distribution (SED). Flux-calibrated lUE spec- 
tra as well as UBVJHK photometry are used to constraint 
log T^ and the amount of extinction. The distances we 

® i-i 

adopted are listed in Table |2] For the Galactic objects, they 

rely on the correlation between the position of the target stars 
and OB associations. Due to uncertainties in these correla- 
tions and in the distance of Galactic associations, the errors 
on the luminosity of the Galactic WNh stars might be under- 
estimated. 

Mass loss rate: it is derived from the absolute strength of 
the emission lines. Ha, Hy8, Hy and He n 4686 are the main 
indicators. The UV resonance lines (Nv 1240, Civ 1550, 
Hen 1640, Niv 1720) are also used. Generally the optical 
and UV wind-sensitive lines give consistent results. 
Terminal velocity: The blueward extent of the UV resonance 
lines is the main indicator. The width of optical lines is the 
secondary diagnostic. 

Surface abundances: The He/H ratio is constrained from the 
relative strength of H and Hell lines. This is different from 
the determination of the mass loss rate in which the absolute 
level of emission is used. The blends He ii/Hy and He ii/Hy6 
together with Hen 4542 and Hen 5412 are the main indi- 
cators. Carbon abundances are determined from C iv 5802- 
5812. Nitrogen abundances result from the analysis of Niv 
4058, Nv 4519-23, Nv 4944, Niv 5205, Niv 5735 and, 
when available, Niv 7103. 

Clumping factor: we used mainly the red wing of the strong 
optical emission lines to determine the volume fillin g factor 
/ (in t he micro-clumping formalism). According to iHillierl 
(Il99ll) . the electron scattering wing depends on the degree 
of inhomogeneity, being stronger in homogeneous winds. In 
practice. He n 4686 provides the best constraints. 
Slope of the velocity field: we have used a value of 1 .0 for the 
so-called p parameter in our modelling of the Galactic tar- 



gets. It provided a good fit of the Balmer and He n emission 
lines. For the LMC stars, we found that fi-l.Q (respectively 
1 .5) lead to a better shape of the optical emission lines of 
Bat 18 (respectively Bat 63). 

The typical uncertainties have been estimated by varying the 
main parameters around the preferred values. When a clear de- 
viation is observed in the fit, we assume we are seeing the ef- 
fect of a one sigma deviation from the best set of parameters. 
The uncertainties are the following: 3000 K on the temperature; 
0.10-0.15 dex on the luminosity; 0.2 dex on the mass loss rates; 
200 km s 'on the terminal velocity; 30 to 40% on the surface 
abundances. 



4. Results 

4.1. SteWar parameters 

The best fit to the observed spectra of our sample stars are shown 
in Figs. [T]to|6] They are usually of excellent quality. Most lines 
are well reproduced. We note a difficulty to reproduce the N v 
4605-4620 features in Galactic stars. They are usually too weak 
compared to the observed spectrum. The other nitrogen lines are 
well reproduced, especially those used to determine the effec- 
tive temperature. The reason for this discr epancy is not known 
at present. Similar problems were noted bv lBouret et al.l(l2012h . 
The derived parameters are gathered in Table |3] 

One of the key parameters we obtained is the hydrogen con- 
tent of the stars. In Fig. |7] we illustrate the determination of the 
He/H ratio in the case of star WR 10. The models all repro- 
duce correctly the Hell 4686 emission. For He/H=0.3, the hy- 
drogen content is too large and the Balmer lines are too strong. 
Alternatively, for He/H=0.7, the hydrogen content is too small 
and the Balmer emission is too weak. Consequently, a value of 
He/H=0.5 is preferred. This also gives an estimate of the uncer- 
tainty on the He/H determination. We see that all WNh stars have 
a non negligible fraction of hydrogen in their atmosphere, while 
the two comparison WN4 stars are hydrogen free. 

Th e four Galactic stars were analyzed by iHamann et al.l 
(l2006h . We usually derive lower temperatures. The difference is 
of about 10000 K for WRIO and WR128, but amounts to 20000- 
30000 K for WR7 and WR 18. It is difficult to increase signifi- 
cantly the temperature of our models for those objects since the 
Fev forest around 1300-1400 A becomes very weak, contrary 
to what is observed. On the other hand, N v 4605^620 is bet- 
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Figure 1. CMFGEN model (red) compared to the observed spec- Figure 3. CMFGEN model (red) compared to the observed spec- 
trum (black) of star BAT 1 8 . trum (black) of star WR 10. 
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Figure 2. CMFGEN model (red) compared to the observed spec- Figure 4. CMFGEN model (red) compared to the observed spec- 
trum (black) of star BAT 63. trum (black) of star WR 128. 



ter re produced at such high temperatures (see also iBouret et al.l 
I2012h . 

The luminosities we derive for the Galacti c stars are con- 
sistent with the results of iHamann et al.l (l2006h . Our mass loss 
rate determination are also in good agreement, with differences 
smaller than 0.2 dex. JHamann et al. ( 2006 ) did not find any ev- 
idence for the presence of hydrogen in WR 7 and WR 18, con- 
sistent with our determinations. For WR 10 and WR 128, they 



quote a hydrogen mass fraction of 0.25 and 0.16 respectively. 
This is again in good agreement with our results. 

4.2. Evolutionary status 

Fig. [8] shows the HR diagram with the position of the sam- 
ple stars indicated by triangles (LMC) and squares (Galaxy). 
The evolutionary tracks are from Geneva an d include rotation. 
Two sets of tracks are displayed: the ones bv iMevnet & Maedej 
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(a) 



(b) 



Figure 8. HR diagram with the position of the sample WN stars shown by triangles (for LMC targets) and square s (Galactic targets). 
Circle s correspond to the comparison Galactic WN4 stars. On the l eft panel, the evolutionary tracks are from iMevnet & Maeden 
(l2005h . while on the right panel they are from lEkstrom et all (l2012h . The bold part of the tracks corresponds to a hydrogen mass 
fraction larger than 0. 1 . 







x[k] 



A[A] 



Figure 5. CMFGEN model (red) compared to the observed spec- Figure 6. CMFGEN model (red) compared to the observed spec- 
trum (black) of star WR 1 8 . trum (black) of star WR 7 . 



(I2005h (left) and the ones bvlEkstromet al.l (l2012h (right). The 
reference solar metallicity is Z=0.020 in the Meynet & Maeder 
tracks, while it is Z=0.014 in the Ekstrom et al. tracks. We show 
in bold the part of the evolutionary tracks for which the hydrogen 
mass fraction is higher than 0. 1 . This value is lower than the val- 
ues we determine for all WNh stars of our sample. The Galactic 
stars cannot be reproduced by the 2005 tracks. Their position in 
the HR diagram is reached only by hydrogen-free tracks. With 
the more recent tracks, the stars are at the limit of the region 



reached by ~25 M0 track still containing a fraction of hydrogen. 
They can be only marginally represented by those tracks. LMC 
tracks with normal (300 km s"') and fast (550 km s"') rotation 
are used to build the HR diagram shown in Fig.|9] The two LMC 
stars (blue triangles) are clearly too hot to be reproduced by stan- 
dard evolutionary tracks (orange lines) which never come back 
to the blue part of the HR diagram with a significant amount of 
hydrogen at their surface, contrary to what Bat 18 and Bat 63 
show. From those comparisons, one can conclude that the LMC 
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Figure 7. Comparison between the observed spectrum of WR 1 Fijure 9. HR diagram with the evolutionary tracks of lBrott et al.l 

(black line) and CMFGEN models with different He/H ratios. (201 1) with an initial rotational velocity of 300 km s"' (orange) 

The lines are H/3 (upper left). Hell 4686 (upper right) and Hy and 550 km s ' (green). The symbols have the same meaning as 

(lower left). The mass loss rate of the models is adjusted so that in Fig. [8] The tracks all correspond to He/H<2.0 (equivalent of 

Hell 4686 is always correctly reproduced. X(H)>0.1). 



WNh stars, and most likely the Galactic ones too, cannot be ex- 
plained by standard evolutiona ry tracks. This is sim ilar to what 
we obtained for the SMC stars (iMartins et al.ll2009l) . 

In this study, we concluded that quasi-chemically homoge- 
neous evolution could be a way to explain the properties of the 
SMC early-type WNh stars. Such evolution takes place when the 
mixing timescale becomes shorter than the nuclear timescale. 
This is possible if the stars rotate initially very fast and do not 
lose too much angular momentum during their evolution. The 
Geneva tracks used so far exist only for an initial rotational ve- 
locity of or 300 km s '. In Fig. |9] we show the HR d iagram 
built with the evolutionary tracks of iBrott et all (l20IIb with a 
LMC composition. Those tracks have been computed for a vari- 
ety of rotational velocities. We have plotted the tracks with ini- 
tial velocities of 300 and 550 km s"'. The latter tracks evolve 
directly to the blue part of the HR diagram from the main se- 
quence. This is very different compared to the 300 km s"' tracks. 
The positions of the stars is well reproduced by the fast rotat- 
ing tracks for initial mass close to 25 M0. In addition, they 
have X(H)>0.1, consistent with our spectroscopic analysis. We 
can thus conclude that the sample stars (at least those from the 
LMC) are consistent with a quasi -chemically homogeneous evo- 
lution. iBestenlehner et al.l dJOUl) studied the WN5h LMC star 
VFTS682 and reached the same conclusion from its position in 
the HR diagra m (they did not constrain the surface abundances). 

Prott et alj (1201 Ih have shown that they could produce this 
type of blueward evolution only at metallicities below that of the 
LMC. Their fast rotating solar metallicity tracks behave more 
classically, with a redward extension from the main sequence. 
But if we look at Fig.[8j we see that the recent Geneva tracks at 
solar metallicity tend to evolve towards the blue quite rapidly af- 
ter the main sequence. The 2005 Geneva models used Z=0.020 
while the 2012 tracks use the revised solar abundances, corre- 



sponding to Z=0.014. This might favour chemically homoge- 
neous evolution at high rotation. And this leaves room for a pure 
homogeneous evolution for higher initial velocities even at solar 
metallicity. 

In Fig. [To] we show the carbon mass fraction as a function of 
nitro gen mass fraction f o r our target stars. The rotating tr acks 
from lMevnet & Maedel (l2005h and lEkstrometaD (1201 2h are 
overplotted. The WNh stars are located in the part of the dia- 
gram corresponding to the minimum of carbon content and max- 
imum of nitrogen content (for their respective metallicity). This 
is typical of stars showing products of the CN cycle at their sur- 
face. This is a very strong indication that the WNh stars are core 
H-buming objects. Indeed, if the stars evolve quasi chemically 
homogeneously and if they show CN equilibrium surface abun- 
dances, this means that their core CN abundances are also at 
equilibrium. Thus, the stars are still in the core H-burning phase. 
If they were more evolved, the central CN abundances would be 
no longer those of CN equilibrium. 

The properties of the sample earl y-type WNh s tars are simi- 
lar to the SMC WNh stars studied bv JMartins et al.l (12009). They 
are relatively unevolved objects. WR 7 and WR 18, the com- 
parison Galactic stars, are as nitrogen rich as the two Galactic 
WNh stars, but show a significantly larger carbon content. This 
is a convincing evidence that they are in a more advanced state 
of evolution since according to the theoretical prediction, after 
reaching a minimum, the carbon mass fraction starts to increase 
due to the onset of core He-burning. These Galactic WN stars 
certainly follow normal evolution and have experienced a clas- 
sical redward evolution in the HR diagram before coming back 
to the left of the ZAMS, after having burnt all of their hydrogen. 
The surface composition of WN stars is thus a key diagnostic of 
their type of evolution. 
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Table 4. Metallicity estimates 




Galaxy 2005 
Galaxy 2012 

LMC 2005 



Star 


d 
[kpc] 


R(kpc) 
[kpc] 


Zpn 


Zhii 


Zhiib 


Zcp 


WR128 


4.30 


7.03 


1.05 


0.97 


1.07 


1.13 


WR18 


3.30 


8.37 


0.88 


0.85 


0.81 


0.95 


WRIO 


4.60 


11.19 


0.60 


0.64 


0.65 


0.66 


WR7 


4.80 


11.30 


0.52 


0.57 


0.60 


0.65 



10 
X(N) [10-3] 



Figure 10. Carbon mass fraction as a function of nitrogen 
mass fraction for the sample stars. The red squares (blue trian- 
gles) are the Galactic (LMC) objects. The Galactic comparison 
stars are shown by the pu rple circles. Evolutionary tracks from 
liVIevnet & IMaedeij (l2005h at solar (LIVIC) metalli city are shown 
by red (blue) lines. The solar metallicity tracks of lEkstrom et al.l 
(1201 2h are shown by the dark red lines. Solid and dashed lines 
correspond respectively to 25 (30 at LMC metallicity) and 60 
M0 stars. 



5. Discussion 

Our results show that CHE takes place in galaxies with aver- 
age metallicities higher than the SMC. Since metallicity gradi- 
ents exist within galaxies, this does not necessarily mean that 
our sample stars are located in region with large metal content. 
Although we have assumed that the LMC stars have a metal 
content 0.5 times the solar value, we have tried to obtain bet- 
ter estimates for the Galactic objects. We have used the posi- 
tion of the Galactic stars in the disk to estimate their metal- 
licity from known O/H and Fe/H gradients established from 
Planetary Nebulae, HII regions and cepheids. Planetary Nebulae 
and HII regions are two types of objects which are commonly 
used to determine the abundance gradient in the Milky Way 
disk. The oxygen abundance, which is easily observed in such 
emission-line objects, is expected to give the interstellar matter 
abundance at the time of the formation of the stars. Cepheids 
provide the iron con tent as a function of galactocentric radius 
dLemasle et al. 2007). The results are given in Table |4] For HII 
regions (IBalser et al..201 li) . we provide two values: one from the 
azimuthally averaged gradient (Zhii) and one from the relation 
established in the azimuthal range corresponding to the star's 
galactic longitude (Zhiib). F or Planetary nebulae, we use the re- 
lation of iHenrv et al.l (l2010l) . The O/H and Fe/H determinations 
give consistent results. From Table|4] we clearly see that WR 1 28 
has a near solar composition. WR 10 is located in a slightly less 
enriched medium, with Z~0.60-0.65. Together with the LMC 
targets assumed to have a metallicity of 0.5 the solar value, we 
conclude that the WNh stars following quasi-chemicaUy homo- 
geneous evolution are found in the interval 0.5-1.0 Zq. 



Notes. Z represents the metallicity. Zpn is obtained from from Planetary 
Nebulae; Zhii and Zhiib from HII regions (see text for discussion); 
Zcp from cepheids. It is express ed in oxygen abunda nce assuming 
12+log(0/H) = 8.66 for the Sun (lGrevesseet"aLl 120071) for planetary 
nebulae and HII regions, and in iron abundance for cepheids. 



We have not determined the metal content of each star di- 
rectly from its spectrum. The Fe v and Fe vi lines in the UV range 
could be used in principle. However, only the Fe v forest around 
1350-1400 A is not dominated by strong emission lines. With 
only a single ionization state available, the determination of the 
Fe content is difficult since temperature uncertainties (of the or- 
der 3000 K) strongly affect the estimate of metallicity. Thus, we 
think metallicity estimates based on galactic gradients are a bet- 
ter than any direct diagnostic we have. 

Our analysis provides good evidence for the occurrence of 
quasi-chemically homogeneous evolution at solar metallicity. 
The current evolut ionary tracks at solar metallicity do not predict 
su ch an evolu t ion dBrott et aI.ll20TTt lEkstrom et al.ll20I2l) . Fig. 5 
of lBrott et al.l (l20Ill) clearly shows that a blueward evolution for 
fast rotating stars is possible only at metallicities below 0.5 Zq. 
This is usually interpreted as a wind effect. At higher metallicity , 
mass loss rates of massive stars are higher (iMokiem et al.ll2007l) 
and thus the removal of angular momentum is easier. Stars brake 
faster and are never fully mixed, the necessary condition to have 
a homogeneous evolution. However, this interpretation assumes 
that the core and the envelope are strongly coupled, so that the 
reduction of the surface angular momentum by stellar wind is 
mirrored in the core angular momentum. If the coupling is weak 
or at l east moderate, the e ffect of winds i s reduced (see disc us- 
sion in lGeorgvet"ai]|20I2h . The models of lBrott et all (|20 III) in- 
clude the transport of angular momentum due to the presence of 
a magnetic field. This strengthens the core-envelope coupling, at 
least partially explaining that the occurrence of quasi-chemically 
homoge neous evolution is corr elated with metallicity and wind 
strength. [Petrovic et al.l (l2005bl) also studied the angular momen- 
tum evolution of single and binary stars and concluded that the 
presence of magnetic field favoured a slowly rotating core^ 

Th e solar metallicity (Z=0.014) models of Ekstrom et al.l 
(l20I2h seen in Fig. [8] (right) show a clear deviation to the blue 
slightly after the zero-age main sequence and above 40 M0. 
These models do not include magnetic coupling. Qualitatively, 
they indicate that quasi-chemically homogeneous evolution is 
not excluded, since they are somewhat intermediate between 
fully chemically mixed models and models with normal mixing. 
Given the uncertainties in the mass loss rates and the physics of 
angular momentum transport, it is thus not excluded that evo- 
lutionary models can reproduce quasi-chemically homogeneous 
evolution. 

According to the present results, to reproduce the properties 
of the early WNh stars in the Galaxy, the following conditions 
are required: 1) a very efficient mixing of chemical species; 2) 
a fast initial rotation; 3) a moderate or weak angular momentum 
coupling between the stellar core and the envelope. In such a 
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way, the effect of stellar winds are reduced to the outer layers. 
This can explain that our WNh stars do not seem to be extremely 
fast rotators. 

The determination of V sin/ for our objects is very diffi- 
cult since most lines are formed throughout the wind and are 
broadened by its expansion. Polarimetry sometimes indicates 
that Wolf-Rayet stars have a flattened shape probably due to ro- 
tation. S ome of these objects may be identified as precursors of 
LGRBs dVink et aLllJoTU iGrafener et al.ll2012l) . However, even 
without polarimetric information, we can determine orders of 
magnitudes for the projected rotational velocity and clearly state 
that V sin/ values larger than 400 km s 'are safely excluded. 
Above such a value the details of some emission line profiles 
(such as Civ 5802-5812) are smoothed out. The surface pro- 
jected velocities of our sample stars are thus of a few tens to 
a couple of hundreds of km s 'at most. Hence, a moderately 
low surface rotational velocity is not necessarily in contradiction 
with a quasi-chemically homogeneous evolution. 

If the angular momentum coupling between the core and 
the envelope is only moderate , a star evolving h omoge- 
neously may form a collapsar (IWooslev et al.l Il993h and a 
long-soft gamma-ray burst. Recent observations reveal that 
such events are not re stricted to low metallicity environments. 
iLevesque et al.l (1201 Oh report the occurrence of a LGRB in 
a galaxy with 12-Hlog(O/H)-9. 0, larger than the solar value. 
Similarlv. lMannucci et al.l ( 1201 lb show that the mass-luminosity 
relation of LGRB host does not stop at solar metallicity: LGRB 
hosts are simply more massive for a given metallicity. Given the 
effects of stellar winds described above, it is likely that lower Z 
favours CHE because braking is less efficient (even in the case of 
weak coupling, there is still an effect). But this type of evolution 
can also take place in solar metallicity environments. The most 
recent compilation of Wolf-Rayet stars by R Crowther Q con- 
tains 433 objects, among which only seven are of spectral type 
WN3h, WN4h or WN5h. If one assumes that such stars have 
similar properties to our sample stars, one may expect only 1 -2% 
of Wolf-Rayet stars to be massive stars evolving quasi homoge- 
neously. This is most likely an upper limit. Indeed, WN3-5h stars 
being core-H burning objects, their number should be compared 
to the number of normal main sequence stars, which is larger 
than the number of Wol f-Rayet stars. In the LMC, the catalog of 
iBrevsacher et al.l (1 1 9991) lists 10 stars with spectral type WN3-5h 
among 134 objects, corresponding to a fraction of 7% of Wolf- 
Rayet stars being hydrogen-rich early-type WN objects. Even if 
a complete census of the Galactic Wolf-Rayet is still missing, it 
seems that the fraction of WN3-5h is larger in the LMC com- 
pared to the Galaxy. Assuming that all WN3-5h stars experience 
CHE, this is consistent with the theoretical prediction that CHE 
takes place more likely - but not exclusively - at low metallicity 
(lBrottetal.1120111) . 

Fully mixed massive stars should be observed in the O star 
phase when they are close to the zero-age main sequence . There 
are several indications that this is the case. iBouret et al.l (l2003l) 
found that the hot 02III(f*) star MPG 355 in the SMC cluster 
NGC 346 is much younger than the bulk of O stars (1 versus 3 
Myr). They argue that it evolves homogeneously, which is con- 
sistent with its relatively high nitrogen enhancement triggered 
by fast rotation. Consequently, they follow isochrones much 
more "vertical" in the HR diagram and their age cannot be de- 
riv ed by standard isochrones. Similar conclusions are reached 
bv IWalborn et al.l (|2004|) concerning ON2III stars. They are lo- 
cated close to the ZAMS in the HR diagram and feature a large 
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N/C ratio. iMokiem et"an (l2007l) studied a sample of LMC stars 
and reported a correlation between high helium enrichment and 
presence of a mass discrepancy (difference between masses de- 
rived from spectroscopy and from the evolutionary tracks). A 
natural explanation is that the use of normally rotating tracks 
for stars actually following homogeneous evolution will lead to 
a strong overestimate of the evolutionary mass. They thus con- 
clude that some of the LMC O stars could be fully mixed. They 
also note that three of their 02 stars are located blueward of the 
main sequence. The indirect evidence for quasi-chemically ho- 
mogeneous evolution in the O star phase are rather compelling, 
at least in the Magellanic Clouds. Similar results for Galactic O 
stars do not (yet) exist. 

An alternative explanation for the position of the WNh stars 
in the HR diagram is binarity. Evolutionary calculations of mas- 
sive binary systems indicate that as soon as the two compo- 
nents interact through mass transfer, the evoluti onary tracks are 
sever e ly modified compared to single star tracks (IWellstein et al.l 
1200 ll: iPetrovic et al.ll2005al) : the accretor can turn back to the 
blue part of the HRD. Recent surveys of young star clusters in- 
dicate that a significant fraction - ranging from 30 to 70 % - of 
massive stars have a companion dSana et al.ll20I2l[20I3l) . Mass 
transfer (and even merging) can oc cur during the early p hases of 
evolution in short-period systems (Ide Mink et al.ll2012b . In such 
a case, the accreting star receives mass and angular momentum. 
It can appear bluer than the bulk of the surrounding population, 
just as blue stragglers in globular clusters. This could explain the 
position of our sample WNh stars in the HR diagram. We have 
seen that the surface composition of the WNh stars was typical of 
CN equilibrium. In the binary scenario, the surface appearance 
of the accreting star depends strongly on the phase at which mass 
transfer occurs. If it is during the expansion of the primary, at the 
end of central hydrogen burning, the material dropped onto the 
secondary can show the signatures of CN processing. But de- 
pending on the initial period of the binary system, mass transfer 
can occur either before or much after hydrogen core exhaustion. 
In that case, no special surface abundance pattern is expected. 
It appears unlikely that binarity can explai n the properties of al l 
the early-type WNh stars we have studied. iFoellmi et al.l (|2003|) 
studied the binary status of Bat 18 and Bat 63 and found no ev- 
idence for radial velocity nor photometric variations, consistent 
with them being single stars. No clear evidence for binarity ex- 
ists for WR 10 and WR 128 either. We thus conclude that bina- 
rity could in principle explain the properties of some early -type 
WNh stars, but the ones we studied here are better accounted for 
by CHE. 

All in all, there are now several robust indications that quasi- 
chemically homogeneous evolution happens at different metal- 
licities, up to solar. An understanding of the conditions under 
which it appears is necessary since the current evolutionary mod- 
els do not predict it at high metallicity. Beyond the understand- 
ing of a somewhat limited number of objects, their study is of 
importance to better constrain the mixing mechanism in massive 
stars. With better predictions, it will be easier to interpret data 
on the occurrence of LGRBs in various environments. 

6. Conclusion 

We have presented the analysis of early-type H-rich WN stars in 
the LMC and the Galaxy. The stellar and wind parameters have 
been determined through spectroscopic analysis in the UV and 
optical range. Atmosphere models and synthetic spectra com- 
puted with the code CMFGEN have been used. The main results 
are: 
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• the early- type WNh stars are hot objects located on the blue 
part of the zero age main sequence. 

• in spite of their Wolf-Rayet classification, they still contain 
a significant amount of hydrogen (X(H) > 0.2). In addition, 
their C and N composition is typical of the CN equilibrium. 
Altogether, this indicates that the WN3-5h stars are core- 
hydrogen burning objects. Comparison to H-free WN4 stars 
indicates clear differences in the evolutionary state of both 
types of objects (WNh versus WN). 

• classical evolutionary tracks cannot reproduce both the po- 
sition in the HR diagram and the chemical composition of 
these objects. Fully mixed models evolving blueward from 
the zero-age main sequence are consistent with such proper- 
ties. WN3-5h stars most likely follow quasi-chemically ho- 
mogeneous evolution 

• estimates of the metallicity in the direct environment of the 
target stars reveals 0.5 < Z < 1.0, indicating that quasi- 
chemically homogeneous evolution happens at solar metal- 
licity. This favours a relatively moderate coupling between 
the core and the envelope for angular momentum transport, 
since the effect of stronger stellar winds at higher metallicity 
does not prevent quasi-chemically homogeneous evolution, 
contrary to the cuiTent predictions of evolutionary models . 
This is also consistent with the "not extreme" rotational ve- 
locity of the WN3-5h stars. 

• the finding that quasi-chemically homogeneous evolution 
takes place in solar metallicity environments and the discov- 
ery of LGRBs in metal-rich galaxies support the scenario 
according to which the progenitors of LGRBs result from 
quasi-chemically homogeneous evolution. 

The total number of WN3-5h stars in the Galaxy is relatively 
small (1-2% of all Wolf-Rayet stars). If all WN3-5h stars fol- 
low quasi-chemically homogeneous evolution, this shows that 
this type of evolution is rare. It is however crucial to understand 
it, not only in the context of LGBR production, but also to bet- 
ter constrain the physical processes in stellar interior, especially 
regarding mixing and angular momentum transport. Future the- 
oretical studies should help improve our knowledge of such ef- 
fects. 
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Table 2. Photometry and distance of the target stars. 



Star 


ST 


U 


B 


V 


J 


H 


K 


distance 
[kpc] 


Galaxy 


WR7 
WRIO 
WR18 
WR128 


WN4 
WN5h 
WN4 
WN4(h) 


11.21 

10.62 

11.10 

9.59 


11.68 
11.30 
11.40 
10.54 


11.40 
10.90 
10.60 
10.46 


9.97 

10.05 

8.57 

9.97 


9.67 
9.89 
8.21 
9.84 


9.27 
9.61 
7.68 
9.62 


3.67 
4.68 
2.20 
9.37 


LMC 


Bat 18 
Bat 63 


WN3h 
WN4ha 


13.57 
13.62 


14.45 
14.45 


14.67 
14.60 


14.63 
14.84 


14.49 
14.62 


14.40 
14.50 


50.12 
50.12 



Notes. References for the distances are: iMassev et al.l (l200lh for WRIO. Ivan der Huchtj (1200 ih for WR18, 7, 128 . The distance modulus to the 
LMC is assumed to be 18.5. 

Table 3. Parameters of the target stars determined from spectroscopic analysis. 



Star 


ST 


[kK] 


T. 
[kK] 


'°s4 


R. 

[Rq] 


log(M) 


[kms-1] 


f 


X(H) 


He/H 


X(C) 
[xlO-^^] 


X(N) 
[xlO-3] 


Galaxy 


WR7 
WRIO 
WR18 
WR128 


WN4 
WN5h 
WN4 
WN4(h) 


60.0 

53.5 
56.0 
57.0 


80.8 
55.2 
74.1 
59.9 


5.40 
5.45 
5.30 
5.50 


2.57 
5.79 
2.73 
5.43 


-4.80 

-5.40/-5.45 

-4.60 

-5.30 


1600 
1400 
2200 
1800 


0.1 
0.1 
0.3 
0.1 


<0.01 
0.33+0.11 

<0.01 
0.26+0.07 


>10 
0.50+0.20 

>10 
0.70+0.20 


14.7+5.0 
7.9+3.0 

17.7+5.0 
6.2+1.6 


11.0+5.0 
14.0+5.0 
10.3+4.0 
11.0+3.0 


LMC 


Bat 18 
Bat 63 


WN3h 
WN4ha 


60.0 

58.5 


72.8 
68.9 


5.50 
5.45 


3.54 
3.73 


-5.02 
-5.45 


1800 
2000 


0.3 
0.1 


0.25+0.10 
0.40+0.20 


0.80+0.20 
0.35+0.15 


5.7+2.0 
20.0+15.0 


6.6+2.1 
2.3+1.6 
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